Strong electron correlation effects in non-volatile electronic memory devices 
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We investigate hysteresis effects in a model for non-volatile memory devices. Two mechanisms 
are found to produce hysteresis effects qualitatively similar to those often experimentally observed 
in heterostructures of transition metal oxides. One of them is a novel switching effect based on a 
metal-insulator transition due to strong electron correlations at the dielectric/metal interface. The 
observed resistance switching phenomenon could be the experimental realisation of a novel type of 
strongly correlated electron device. 
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The inclusion of strong correlation effects in semi- 
conductor electronic devices has long been a goal of 
solid state physicists. Significant progress has been 
achieved in recent years in studies of metal-insulator- 
metal (MIM) structures that display different types of re- 
sistance switching by the application of a voltage or cur- 
rent pulse; a prototype of two-terminal nonvolatile ran- 
dom access memories. The insulator of the MIM struc- 
ture is either a transition metal oxide(ll,t^,^,^,^] or an 
organic materially, LDli in which the electron correlations 
usually play an important role. The switching is directly 
related to the large hysteresis observed in the J- V charac- 
teristics. These systems are collectively called resistance 
random access memories (RRAM) and the most funda- 
mental current issues for this technology are the differ- 
ent hysteretic behaviour observed in the current-voltage 
characteristics and their relation to non-volatility. 

The goal of the present work is to investigate the qual- 
itatively different hystereses of I-V curves in the MIM 
structures that have been experimentally observed so far, 
and to associate each hysteresis type to an appropriate 
physical mechanism. Our starting point is a recently in- 
troduced model for the RRAM devicefSj. We shall first 
study the basic model predictions for the I-V characteris- 
tics, and then show how it can be extended, incorporating 
the physics of strongly correlated effects. 

The model assumes the existence of an insulating (and 
inert) medium with a nonpercolative structure of metallic 
domains that might correspond to defects, grains, phase 
separated regions, and so onQ. Smaller domains are 
closer to the electrodes and are called "top" and "bot- 
tom" . A large domain occupies the bulk of the system 
that is called "central". Carriers tunnel between do- 
mains under the action of an external electric field (or 
applied voltage) . The probability of charge transfer de- 
pends on phenomenological parameters such as the tun- 
neling rates, the number of states in the domains and 



their occupation (see Ref. |3| for details). The model is 
mathematically defined by the following system of rate 
equations, 

-^^TfN^n^l ~ n>^)U{V) - r^r.,^A^,(l - n^)fUV){l) 

i 

flTI 

-±^TfN,n^{l - n'^UV) - TTn\N,{l - n^)fte{V){3) 

where F" denote the tunneling rates between the elec- 
trodes and domains, and a,/3 = e,t,c,b denote "elec- 
trode" , "top" , "central" and "bottom" . Na is the total 
number of states in electrode or domain a, and n" is 
the occupation. fai3{V) describes the dependence of the 
transfer probabilities between a and /? on the given ex- 
ternal voltage protocol V{t). The subindex i labels the 
domains. 

To gain a qualitative understanding of the model, we 
make some simplifying assumptions. We ignore disorder 
effects and take all "top" and "bottom" domains to be 
identical. The functions fafsiV) are assumed to be in- 
dependent of a and /S. We also take a unique tunneling 
rate for the two interfaces F*^* — F'"^ = F*^^* and similarly 
inside the dielectric F*'^ = F'^'' = F™*. Therefore, we 
can just consider the behaviour of the "average" top and 
bottom domain occupations, that we denote n* and n''. 
Finally the metallic electrodes are assumed half-filled. 

The set of equations simplifies to three coupled non- 
linear differential equations. 



~dt' 
dn'^ 

~dt 
dn* 



- [F^-*^(l_n^)-r"*n^A^,(l-n^)]/(T/) (4) 
= [T"''Nbn\\ ~ n"") - V'n''Nt{l - n')]f(V) (5) 
= [F™*iVcn^(l-n*)-F'=^*n*— ]/(y) (6) 



that have to be solved for the unknown n"(i) with 
a = t,c, b. The coefficients f{V) are time dependent 
through the appUed voltage protocol V{t). There is no 
general solution for this model, but specific cases can 
be analysed. As the central domain is taken to be large 
compared to top and bottom, we may assume that its oc- 
cupation remains approximately constant, i.e., -4- « 



and n'^(t) 



Under this condition, the system of 



Eqs. (4)-(6) can be formally solved and the result for the 
time dependent bottom domain occupation is. 



n\t) 



-^-IoP(t')dt' 



Q{t')eloPi'")'^*"dt' + n\0)) (7) 



where P{t) = /(y)[r™*iVc(l - <) + T'^^'iVe/S] and 
Q{t) = /(T/)[r^^* Are/2]. The occupation of top is simply 
given by, n*{t) = 1 — [(1 — n^) / n'^]n'' (t) . Taking the time 
derivative on both sides we see that n^ controls the charge 
accumulation in the model, that is, [dn*-/dt — dn^/dt]. 

Equation Q provides further insight into the non- 
volatility properties. As the system is strictly always out 
of equilibrium under applied voltage, the key quantity 
to consider is the relaxation time under either applied 
bias Vb or during a voltage pulse Vp {Vb <^ Vp). The 
main dependence of the relaxation time comes from the 
first exponential factor of (0. With the rather natural 
assumption that the tunneling probability dependence is 
f{V) ~ e^, one finds that the ratio of relaxation times 
is exponential in (Vp — Vb). Thus bias relaxation can 
be made orders of magnitude slower than the switching 
time, which is consistent with the observed behavior in 
actual systems. We further note that when the device is 
disconnected from the external battery there remains an 
internal voltage Vint due to the inhomogeneous charge 
distribution. The coulombic forces would produce the 
subsequent relaxation of this state, and this process is 
crucial for the long time non-volatility characteristics of 
a system Ig. The previous analysis implies the existence 
of a compromise: driving the system to a strong charge 
inhomogeneous state (with a strong or long Vp pulse) 
would provide a desired large set/reset resistance ratio 
but, on the other hand, would also lead to faster relax- 
ation, i.e., shorter non- volatility time. These results are 
quite consistent with the investigations of relaxation in 
Nb-doped SrTiOa/LaTiOs.s/Au heterostructures0. 

Let us now turn to the predictions of the model for 
the I-V characteristics and hysteresis effects. The ap- 
plied voltage protocol is ^ Vmax — > —Vmax — * with a 
constant sweep velocity. For simplicity, all three domains 
are assumed to be initially half filled. Following Simmons 
and Verderberjl^, we take the function f(y) — sinh(fcy) 
where k depends on various material parameters, thus we 
set k — 1. The (positive) current in and out of the system 
is given by the carriers entering the bottom domain and 
leaving the top domain respectively. In Fig.^ we show 
the results for two qualitatively different hystereses that 




FIG. 1: Current-voltage characteristics. Nt = Nb — lO'', 
iVc = 10^° and F^"* = 2xl0"^ F-^^* = 6x10"" (left panel) or 
pint ^ 8xl0-^ F''^* = 6x10"^^ (right panel). The V protocol 
cycle begins at the thick arrow. 



the model predicts depending on the choice of parame- 
ters. Their main difference is that one shows low to high 
current switching (left) while the other shows the oppo- 
site behavior (right). It is not difficult to understand the 
reason for this contrast. The current out (for current in, 
a similar analysis can be done) depends mainly on the 
filling level of the top domains. The larger their occupa- 
tion, the more carriers that are available to tunnel out. 
During the positive side of the voltage protocol, V is ini- 
tially increased and then decreased down to zero. During 
this phase, charge is transferred from the center domain 
to the top and from there to the electrode. Therefore, 
if more charge enters the top domain than leaves it, the 
filling level during the increasing positive-^ ramp would 
be lower than the filling level during the subsequent de- 
creasing positive-^ ramp. This would lead to low-to- 
high current switching as shown in the left panel. There- 
fore, the qualitative form of the hysteresis depends on 
the intrasystem (top/center and bottom/center) charge 
transfers compared to the interface (top/electrode and 
bottom/electrode) ones, both of which are controlled 
by the model parameters (basically the FA^ products). 
Assuming larger interface transfers should be appropri- 
ate for systems such as Au/SiO/Au studied in the pio- 
neering work of Simmons et a-/.|10l lll|. where hysteresis 
and negative differential resistance (NDR) were observed 
and qualitative agreement with our model is found, as 
shown in the right panel of Fig.^ On the other hand, 
larger intrasystem transfer is a likely assumption for the 
heterostructures with transition metal oxide dielectric, 
where electronic inhomogeneous states are often realized. 
It would also be the case for the system of artificially cre- 
ated domains of Al-nanoclusters, where, in fact, hystere- 
sis similar to the results shown in the left panel of Fig.Q] 
has been observed [TJ. 

While the present model provides useful insight for 
the analysis of relaxation times and basic switching 
mechanisms jg|, the qualitative predictions for hystere- 
sis effect illustrated above show a notable limitation; 
namely, the I-V curves never cross at the origin. 
Experimentally, crossing I-V characteristics are often 
observed[lJ,y, yllsllg, thus providing an important mo- 




FIG. 2: Current- voltage characteristics. Nt = Ni, — 10®, 
N, = 10^0, F"* = 6xl0-^ r""' = 10"" and n*(0) = n^O) = 
n°(0) = 0.1. Basic model (top left), model with interface 
charge dependent tunneling (top left), idem with Schottky ef- 
fect (bottom right), model with occupation driven Mott tran- 
sition and A/T — 3 (bottom right). The V protocol cycle 
begins at the thick arrow. 



tivation to consider two extensions of the basic model, 
Eqs. (l)-(3), that we shall describe next. 

The first one is to incorporate the physics of inter- 
face tunneling rates that depend on the charge accumu- 
lation, i.e., the dependence of V^^^ with n* and n^. The 
underlying physical mechanism is similar to the depen- 
dence of the depletion width with impurity concentration 
in ohmic contacts of metal/semiconductor interfaces |l2j. 
Basically, the accumulated charge density at the interface 
produces additional band bending that leads to a reduc- 
tion in tunneling, thus enhancing r*^^*. Borrowing from 
standard junction tunneling theorvjl2j we can qualita- 
tively model this dependence as Fj'l* = e"'^/^"*' +ir'^^* 
(Ref. 13]) so the tunneling is enhanced with increasing 
charge accumulation in the small interface domains right 
up to the saturation value F*^^* for fully occupied do- 
mains, i.e., n*'''=l. This mechanism is widely assumed to 
be relevant for many resistance switching phenomenona, 
however it is also usually linked to poor non-volatility 
characteristics. 

We incorporated this effect in our basic set of Eqs. (4)- 
(6) and solved for the model behavior. The results are 
shown in Fig. [3 (top right) along with the predictions of 
the original basic model (top left) for the same parameter 
set (for reference). We see that, in fact, this additional 
assumption is sufficient to turn the hysteresis into the 
crossing type, a result that can be understood through a 
similar analysis to that we described earlier|^. An im- 
portant point to mention here is that the hysteresis effect 
does not necessarily need to take place at both electrodes 
simultaneously, since the large center domain acts as a 
buffer that effectively decouples the two interfaces. Hence 
this mechanism would also be relevant for systems with 
two different interfaces where the switching properties 



are controlled by one of them, e.g. in heterostructures 
such as Au/PbTiOs/Lao.sSro.sCoOa where the Au inter- 
face forms a Schottky barrier but Lao.sSro.sCoOs does 
notpj. To investigate a possible effect of Schottky bar- 
riers in our results, we modified the voltage dependent 
part of the tunneling function to the form f{V) — e^-1 
(forward bias) that produces Schottky type I-V charac- 
teristics. The results are shown in Fig.El (bottom left) 
and we find no qualitative differences in the hysteresis. 

As we mentioned before, this mechanism is often 
linked with poor non- volatile characteristics, possibly be- 
cause it lacks a stabilisation effect that could lengthen 
the lifetimes of the charge-accumulated (or -depleted) 
metastable states. A physical mechanism that would en- 
able such stabilisation is likely to be due to a coherent 
effect in which many particles act in a correlated man- 
ner. Here we propose that such a mechanism is based 
on a Mott metal-insulator transition taking place in the 
domains when they become close to half filling. In fact, 
this occupation (or doping) driven metal-insulator tran- 
sition is a strong correlation effect found throughout the 
transition metal oxide series|l4J and even in some organic 
materials [l5|. For concreteness, we shall assume that the 
domains undergo a Mott transition when they are close 
to half-filling. On becoming insulating, the domains open 
an insulating gap A in their excitation spectrum which 
we model by assuming a reduction in the density of states 
at the Fermi energy given by D = e^^^^Do, T being the 
temperature. The dependence of the density of states at 
the Fermi energy was entered implicitly in our equations 
through a multiplicative factor (set to unity) in the func- 
tion f{V), since we were always concerned with tunneling 
between metallic domains. We now consider this effect 
explicitly by modifying f{V) in the model Eqs. (4)-(6) 
to fmit{V) — e^^/^/(F) if the occupation of the cor- 
responding domain is within 10% of half-filling, or else 
fmit{V) = fiV). The results are shown in Fig. El (bot- 
tom right). The results can be understood as follows. 
Initially, the domains are assumed to be lightly filled, i.e. 
in the metallic state. As the applied voltage is increased, 
charge is accumulated at the top domains thus increasing 
their occupation level. Eventually, when they approach 
half-filling, they become insulating and open a Mott gap. 
The tunneling probability to and from the top domain 
gets dramatically reduced, with a consequent drop in the 
output current and also in the rate of further charge in- 
jection. With the reduced charge injection, the occupa- 
tion level of the domain remains about constant and the 
current is low for the rest of the positive voltage sweep 
{-^ Vmax -^ 0). This is also the case for the beginning of 
the negative voltage sweep, as the leakage of charge out of 
the top domain is also greatly suppressed [thus providing 
an enhanced non-volatility of the high resistance state). 
Eventually a large negative V drives enough charge out 
of the top domains and once again they become metallic, 
their correlation gap closes and the system switches back 













1 1 1 1 1 1 1 




40 




' ' 


' "]"•-- 


2e^05 


3D 






/ 1 


- 




t 




^ 


^^ J 


- 




JZ^^^ 


- — 


le+05 


-20 


y. 


^ 




_ 




-4Q 






, 


yj : 


















3 0.5 1.0 1.5_ 


2^;:^--^ 


-50000 


^ 


" 


« — 


-le+OS 


- ''/ 


//- 




- 


,5e+05 


- / 






- 


-2e+05 








- 


..5tH-05 










5 


4 


3 -2 


-1 1 


2 3 4 5 



FIG. 3: Current-voltage characteristics for the model with 
occupation driven Mott transition. The model parameters 
are Nt=Nt=10'^, Nc=10^ and p"*=10-", r'^^^^SxlO"^". 
The current is in carriers per unit of time (ie, Monte 
Carlo step). Inset: experimental I-V characteristics in 
Au/SrTiOa/SrRuOs from Ref. [H. 




10%), consistent with the fact that the systems are nom- 
inally empty band insulators with only light doping due 
to the Cr substitution. Moreover, the intriguing sudden 
high to low current switching at positive applied bias in 
the experiment is a feature that naturally emerges from 
our domain-doping driven metal to insulator transition 
scenario. In Fig. 0] we make more detailed comparison to 
the experimental data by plotting the I-V characteristics 
in different scale types, all of which show excellent qual- 
itative agreement with experiment. We note that while 
the current is proportional to V at low voltages the sys- 
tem is never truly ohmic in the sense that there is no 
actual sample current, rather a series of very incoherent 
charge transfer processes. In fact, current conduction in 
these systems is known to be non-filamentaryJ2|. 

To conclude, we have proposed a model for resistance 
memory switching and find the different physical mech- 
anisms that lead to various qualitatively different hys- 
teresis effects in the I-V characteristics. We proposed 
a new scenario for memory switching in strongly corre- 
lated transition metal oxide heterostructures, based on 
a metal-insulator transition driven by charge injection 
(doping) to domain structures in physical proximity to in- 
terfaces. The information (high resistance state) is stored 
by driving the systems to a qualitatively different state 
which, at the same time, provides a novel mechanism for 
the enhancement of the non- volatility. Our model results 
are supported by excellent qualitative agreement with ex- 
perimental data in various heterostructures. 
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FIG. 4: Idem as in Fig.|3 in normal scale (top), semi-log 
scale (middle) and log-log scale (bottom). The left panels 
show experiments in Au/SrZrOs/SrRuOs from Ref. U- 



to a high current state. 

In order to give further support to this mechanism we 
performed Monte Carlo simulations |3 on the model given 
by Eqs. (l)-(3) supplemented with the Mott transition 
effect. We looked for parameters that qualitatively re- 
produce the experimental results in the MIM structures 
with Cr-doped SrTiOa or SrZrOs of the IBM group[lE3- 
Our results are shown in Fig. Inland the experiments are 
shown in the inset. (We modified the voltage protocol 
to -Vmax — > Vrnax "^ -Vmax to match the experiment.) 
Not only is the qualitative agreement very satisfactory, 
several details are also worth pointing out. The initial oc- 
cupation level of the domains that allow best comparison 
to the experimental data was found to be small (around 
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